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Thermal  management  is  crucial  for  improving  the  charge-discharge  efficiency  and  cycling  life  of  lithium 
ion  battery.  In  this  paper,  a  mathematical  model  coupling  electronic  conduction,  mass  transfer,  energy 
balance  and  electrochemical  mechanism  is  developed.  Lithium  ion  diffusivity  and  chemical  reaction  rate 
of  cathode  material  are  estimated  by  comparing  simulated  results  with  experimental  data  of  pulse  test  at 
various  current  charge-discharge  rates  (0.2C,  0.5C,  1 C,  2C)  and  operating  temperatures  (0  °C,  1 0  °C,  25  °C, 
55  °C).  The  modeling  results  are  further  validated  in  aspects  of  electrochemical  performance,  thermal 
performance  and  electrochemical-thermal  coupling  effects,  which  show  well  agreement  between  the 
modeling  results  and  experimental  results.  The  modeling  results  show  that  lithium  ion  concentration 
gradient  in  both  liquid  phase  and  solid  phase  are  greatly  affected  by  temperature,  and  the  lithium  ion 
concentration  gradient  increase  when  temperature  decrease.  This  phenomenon  results  in  the  capacity 
losses  and  power  loses  of  lithium  ion  battery  during  low  temperature  operation.  The  reversible  heat 
generation  during  charging  process  is  equal  with  the  heat  consumption  during  discharging  process.  It  is 
also  indicated  that  the  reversible  heat  is  dominant  at  low  rate  discharging  process  and  irreversible  heat 
is  dominant  at  high  rate  discharging  process.  Proper  cooling  system  should  be  added  to  keep  battery 
temperature  within  safety  range  during  high  rate  current  charging/discharging. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium  rechargeable  batteries  of  various  types  have  been  mov¬ 
ing  rapidly  toward  commercialization  ascribing  to  their  potential 
advantages  in  power  density,  cost,  safety,  performance,  or  design 
flexibility.  Recently,  it  has  been  paid  more  and  more  attention  to 
the  application  on  electric  vehicle  (EV)  and  hybrid  electric  vehi¬ 
cle  (HEV).  The  designers  need  a  complete  understanding  of  the 
characteristics  of  lithium  ion  battery  [1  ].  However,  choosing  among 
the  various  chemistries  and  optimizing  the  design  of  battery  for  a 
certain  application  field  will  be  a  costly  and  time-consuming  exper¬ 
imental  undertaking.  Numerical  modeling  and  simulation  of  the 
battery  is  crucial  for  improving  the  fundamental  understanding  of 
battery  operation  and  for  designing  appropriate  thermal  manage¬ 
ment  system. 

Several  electro-thermal  models  of  lithium  ion  battery  have 
been  presented  in  published  literature.  Newman  and  Pals  [2,3] 
considered  isothermal  and  adiabatic  discharge  behaviors  in  the 
thermal  models  of  cell  and  battery  stack  by  adding  an  energy  bal¬ 
ance  equation  to  the  pseudo  two  dimensional  porous  electrode 
models  of  Doyle  et  al.  [4].  Other  researchers  [1,5,6]  further  devel¬ 
oped  electrochemical-thermal  model  by  adding  separated  heat 
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source  terms  and  by  employing  temperature  dependent  param¬ 
eters  (lithium  ion  diffusivity,  conductivity,  etc.).  Generally,  heat 
generation  within  lithium  ion  battery  can  be  classified  into  three 
types:  (1)  the  activation  polarization  irreversible  heat  due  to  the 
electrochemical  reaction  polarization  between  active  material  par¬ 
ticle  surface  and  the  electrolyte;  (2)  ohmic  heat  due  to  the  ohmic 
potential  drop  and  (3)  reversible  reaction  heat  due  to  entropy 
change  during  charge-discharge. 

Furthermore,  for  a  given  mathematic  model,  although  most  of 
the  model  parameters  are  determined  from  experiments  or  from 
the  literature,  others  are  only  estimates.  Some  of  the  parame¬ 
ters  are  adjusted  to  ensure  good  agreement  between  the  model 
results  and  the  experimental  data,  for  example  kinetic  and  trans¬ 
port  properties.  The  uncertainty  in  the  estimated  parameter  will 
surely  influence  the  model  reliability  and  accuracy.  The  more  the 
parameters  are  determined  independently  through  experiments  or 
characterization  techniques,  the  more  robust  the  model.  In  addi¬ 
tion,  the  comparison  between  experimental  charge  and  discharge 
data  and  modeling  results  helps  battery  designer  deeply  under¬ 
stand  the  effects  of  various  parameters  (thermodynamic,  kinetic, 
and  design)  on  the  battery  performance  with  different  operating 
conditions,  such  as  charge-discharge  rate  and  temperature.  Thus, 
the  modeling  efforts  can  possibly  reduce  the  experimental  work 
required. 

In  this  study,  an  electro-thermal  coupling  model  of  lithium 
ion  battery  is  proposed.  Model  parameters  which  are  sensitive  to 


228 


Y.  Ye  et  al.  /  Journal  of  Power  Sources  199  (2012)  227-238 


Nomenclature 

Acell  area  of  the  positive  (both  sides)  that  has  opposing 
negative  (m2) 

Ci  f  reduced-lithium  in  active  material  ( mol  nrr 3 ) 

Ci  <max<i  maximum  concentration  ( mol  nrr 3 ) 

^i, surf, i  lithium  ion  concentration  on  active  material  parti¬ 

cles  (mol  nr3) 

CPti  heat  capacity  (J  (kg  K)"1 ) 

Di  j  solid  phase  diffusivity  (m2  s_1 ) 

Dio,/  solid  phase  diffusivity  at  reference  temperature 

(m2  s_1) 

Ea,D,i  active  energy  for  diffusion  (kj  mol-1 ) 

Ea  k  i  reaction  active  energy  (kj  mol-1 ) 
h  heat  transfer  coefficient  (W  nr2  K-1 ) 

I  cell  current  (A) 

/app  cell  current  density  based  on  Acell  (A  nr 2 ) 

jo, /  exchange  current  density  (A  nr 2 ) 

Jioc,/  local  current  density  (A  nr 2 ) 

/<0,2  reaction  rate  constant  (m2-5 mol-0,5 s-1 ) 

kf  thermal  conductivity  (W  (m  I<)-1 ) 

I/  thickness  (f±m) 

Qact  active  heat  generation  (J  nr 3 ) 

Qohm  ohmic  heat  generation  (J  nr 3 ) 

Qrea  reaction  heat  generation  (J  nr 3 ) 

r  radius  distance  variable  of  particle  (m) 

Kceii  total  resistance  in  battery  (Q) 

Ri  characteristic  radius  of  electrode  particles  (pan) 

Sa  (  specific  surface  area  (nr1 ) 

SOCq  i  initial  state  of  charge 

t  time  (s) 

t+  transference  number  of  Li  ion  species  dissolved  in 

liquid 

T  absolute  temperature  (K) 

^amb  ambient  temperature  (I<) 

Uf  thermodynamic,  open  circuit  voltage  (V) 
v  thermodynamic  factor  relating  to  electrolyte  activ¬ 

ity 

x  distance  variable  through  a  cell  component  (m) 

y  dimensionless  radial  distance  of  particles 

x  average  composition  variable  of  Li  in  Li^C6 

y  average  composition  variable  of  Li  in  LiyMn204 

Greek  letters 

aai  transfer  coefficient  for  anodic  current 
otc j  transfer  coefficient  for  cathodic  current 

£1(I-  active  material  volume  fraction 

£2,2  volume  fraction 

02  electric  potential  (V) 

Yi  Bruggeman  tortuosity  exponent 

k  ionic  or  electronic  conductivity  (S  nrr1 ) 

Pi  density  (kg nr3) 

o'/  solid  phase  conductivity  (S  m- 1 ) 

Subscripts,  superscripts  and  acronyms 
0  initial  or  equilibrated  state 

1  solid  phase 

2  liquid  phase 

amb  ambient  (temperature) 

app  current  applied  to  battery 

Cell  battery 

Ch  charge 

Disch  discharge 

irr  irreversible 


n  negative  electrode 

p  positive  electrode 

re  reversible 

ref  reference  composition  of  relative  to  a  Li/lithium  ion 
reference  electrode  or  reference  temperature  for 
Arrhenius  formula 
s  separator 

surf  surface  of  active  material  particles 


temperature  change  are  estimated  by  comparing  modeling  results 
with  experimental  data.  Temperature  sensitive  parameters  (such  as 
diffusivity  of  lithium  ion  in  solid  and  liquid  phase  and  reaction  rate 
constant)  are  estimated.  After  comprehensive  model  validation  at 
different  operating  temperatures  and  charge-discharge  rates,  the 
influence  of  operating  temperature  on  lithium  discharge  behavior 
and  lithium  ion  concentration  distribution  are  simulated  and  dis¬ 
cussed.  The  roles  of  heat  sources/sinks  of  electrochemical  reaction 
in  different  operating  conditions  are  analyzed. 

2.  Experiments 

Commercial  LiMn204  battery  (90  mm  x  57  mm  x  25  mm  pris¬ 
matic  battery,  capacity  11. 5 Ah,  LiC6  anode  material,  LiMn204 
cathode  material,  1 .2  mol  L-1  LiPF6  in  PC/EC/DMC  solvent)  is  exper¬ 
imentally  characterized.  A  LAND  (CT2001-B,  China)  test  system, 
8-channel  A  to  D  converter  and  a  computer  data  logger  is  used 
to  monitor  the  charge-discharge  current  and  battery  voltage.  Real¬ 
time  graph  is  created  and  stored  in  computer.  The  battery  is  placed 
in  a  temperature  control  box  for  being  operated  at  different  tem¬ 
peratures. 

The  full  capacity  of  the  battery  is  measured  at  25  °C  with  follow¬ 
ing  steps:  (1)  the  battery  is  fully  charged  to  3.8  V  with  a  constant 
current  of  2  A  and  with  constant  voltage  at  3.8  V  till  charge  current 
declined  to  200  mA  at  25  °C;  (2)  according  to  the  method  declared 
by  the  manufacturer.  Then,  battery  is  discharged  with  a  constant 
current  of  2  A,  down  to  2.5  V,  the  end-of-discharge  voltage  specified 
by  the  battery  manufacturer. 

The  initial  OCV  of  batteries  was  measured  using  multimeter  with 
5.5  digit  resolution.  The  battery  is  discharged  at  constant  current 
(0.2C,  0.5C,  1C,  2C)  at  room  temperature  (25  °C),  so  as  to  character¬ 
ize  the  electrochemical  behavior  of  lithium  ion  battery  at  normal 
temperature. 

A  temperature  sensor  probe  was  stuck  on  the  center  site  of 
the  batteries  to  monitor  the  surface  temperature  change  during 
operation.  The  batteries  were  carefully  wrapped  by  thermal  insu¬ 
lation  material-cotton.  Another  temperature  sensor  measured  the 
ambient  temperature  in  situ.  Lastly,  the  batteries  were  charged  and 
discharged  at  some  constant  current  (0.2C,  0.5C,  2C)  continuously. 
Between  the  charging  process  and  discharging  process,  10  h  were 
given  to  let  the  batteries  temperature  down  to  the  ambient  tem¬ 
perature.  During  full  process,  battery  surface  temperature  will  be 
measured  and  recorded  in  situ  by  multi-channel  temperature  log¬ 
ger.  We  tried  to  conduct  these  experiments  in  insulating  condition 
so  as  to  eliminate  the  influence  of  ambient  temperature,  however, 
the  experimental  condition  cannot  be  ideally  insulate,  so  it  is  actu¬ 
ally  a  quasi  insulating  condition  with  relatively  small  heat  transfer. 

Pulse-test  at  different  ambient  temperatures  (0°C,  10°C,  25  °C, 
55  °C)  and  different  currents  (0.2C,  0.5C,  1C,  2C)  are  performed  to 
characterize  dynamic  performance  of  lithium  ion  battery.  It  should 
be  noted  that  the  temperature  dependent  parameters,  such  as 
lithium  ion  diffusivity  in  solid  and  liquid  phases  and  reaction  rate 
can  be  extracted  from  these  experimental  data. 
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3.  Model  development 

3.1.  Model  assumptions  and  calculation  domain 


battery,  Eq.  (3);  at  boundary  2  and  boundary  3,  there  is  no  charge 
flux,  Eq.  (4)  and  the  boundary  condition  is  set  as  isolation 


Keff 

I  =0; — Vv0 1 

^  I  x=0  [2  ^ 


y^eff 


=  0; 


x=Ln  +Ls 


x=Ln  +LS  +Lp 
^eff 

-ijv*, 


=  0 


X=Ln 


(3) 

(4) 


3.2.12.  Solution  phase.  The  governing  equation  for  electronic 
charge  balance  in  solution  phase  is  expressed  as 


—  V02  + 


9  In  / 
9  In  c2 


1(1-4) 


Vc2 

c2 


} 


—  Sa,JlocC 


(5) 


Kf  =  K28$ 


(6) 


Liquid-junction  potential  is  introduced  in  Eq.  (5)  with  expres¬ 
sion 


„  2 RT„  31n  / .  N  2RT 


(7} 


where  v  is  the  thermodynamic  factor  relating  to  electrolyte  activity, 
it  is  temperature  and  concentration  dependent. 

There  is  no  flux  at  external  boundaries  (boundary  1  and  bound¬ 
ary  4),  Eq.  (8);  (p2  is  taken  to  be  continuous  at  inner  boundaries 
(boundary  2  and  boundary  3). 


902 

dx 


x=0 


902 

dx 


=  0 


x=Ln  +Ls+Lp 


(3) 


The  electro-thermal  model  for  LiMn204  lithium  ion  battery  in 
this  paper  is  developed  based  on  the  pseudo  two-dimensional 
model  of  Doyle  et  al.  [4,7].  The  main  model  assumptions  were 
shown  as  follows: 

(1)  Gas  generated  during  operation  is  neglected,  and  only  liquid 
phase  and  solid  phase  in  lithium  ion  battery  are  considered. 

(2)  Side  reactions  during  operation  are  neglected. 

(3)  Active  materials  in  solid  electrodes  are  considered  to  be 
homogenous,  and  are  composed  with  spherical  particles. 

(4)  The  effects  of  current  collectors  on  lithium  ion  transfer  and  heat 
transfer  are  neglected. 

Fig.  1  shows  a  schematic  calculation  domain  of  one-dimensional 
(ID)  battery  model  [4,7,8].  Two  inner  boundaries  (anode/separator 
interface  boundaries  2  and  cathode/separator  interface  boundaries 
3)  and  two  external  boundaries  (anode/current  collector  interface 
boundary  1  and  cathode/current  collector  interface  boundary  4)  are 
shown  in  the  figure. 


3.2.2.  Mass  balance 

3.22.1.  Solid  phase.  The  mass  balance  of  lithium  ions  in  an  inter¬ 
calation  particle  of  electrode  active  material  is  described  by  Fields 
law.  Let  the  variable  y  equal  to  r/Rif  the  mass  transport  within  solid 
phase  can  be  described  as 


y2R, 


dci  j  9 

~dt  +  9y 


(9) 


where  Rt  is  the  particle  radius.  The  lithium  concentration  at  the 
surface  of  the  particles,  y  =  1,  is  coupled  to  the  concentration  and 
flux  in  the  1 D  model  for  the  charge  and  material  transport  in  the 
electrolyte.  The  flux  at  the  center  of  sphere,  y  =  0,  is  set  to  zero, 
because  there  is  no  species  source. 


3.2.22.  Solution  phase.  Solution  phase  material  balance  for  LiF6 
dissolved  in  the  liquid  phase 


D? 


Li$a,dloc,i 


(1-4) 


(10) 


3.2.  Governing  equations 


3.2.1.  Electronic  charge  balance 

32.1.1.  Solid  phase.  Electronic  charge  balance  for  solid  phase  can 
be  expressed  as  follows: 


LiV- 


y-eff 


£ 2  ^01  )  —  ^i^a,d\oc,i 


Sa,i=^f--.Kf  =  K  Ifip 


(1) 

(2) 


where  0 1  is  arbitrarily  set  at  zero  at  boundary  1 ;  at  boundary  4,  the 
charge  flux  is  set  to  be  equal  to  the  average  current  density  of  the 


Df  =  D2s %  (11) 

The  flux  of  liquid  species  is  set  to  zero  at  external  bound¬ 
aries  (boundary  1  and  boundary  4),  liquid  species  flux  and  species 
concentration  are  taken  to  be  continuous  at  inner  boundaries 
(boundary  2  and  boundary  3). 

3.2.3.  Electrochemical  kinetics 

The  local  current  per  active  material  area  is  calculated  using  the 
Bulter-Volmer  equation: 
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where  ji0C  I  is  driven  by  overpotential,  pi,  defined  as  the  difference  4.  Model  calibration  and  validation 
between  solid  and  electrolyte  phase  potentials  minus  U ,  the  ther¬ 
modynamic  equilibrium  potential  of  the  solid  phase.  4A.  Model  calibration 


Vi  —  01,  i  02,  i  ^ref  ,i 


(13) 


where  the  equilibrium  potential,  U,  is  taken  to  be  a  function  of 
the  solid  phase  lithium  ion  concentration  at  the  particle  surface.  In 
Eq.  (14),  exchange  current  density,  jo  i,  acts  as  a  bridge  connecting 
concentrations  in  both  solid  phase  and  liquid  phase: 


Jo,;  =  FkiC%°,i  (c1>maXii  -  c1>surf>I-)"OJc"c, 


surf,  i 


(14) 


where  /q  is  a  reaction  rate,  considered  temperature  dependent  in 
this  paper.  aa  and  ac  are  the  anodic  and  cathodic  transfer  coeffi¬ 
cients,  respectively. 


3.2.4.  Energy  balance 

The  following  equation  governs  energy  balance  in  lithium  ion 
battery.  There  are  three  parts  of  heat  sources  during  charge  and 
discharge  processes,  including  reaction  heat  Qrea,  ohmic  heat  Qohm 
and  active  polarization  heat  Qact. 

UPiCp/A  +  V  •  (^Vr)  =  L,(Q rea  +  Qact  +  Qohm)  (15) 

where  active  polarization  heat  is 

Qact  =  Sa,d loc,i(01,i  —  02, i  ~  Ui)  (16) 

Reaction  heat  is 

Orea=Sa,Jloc>IT^  (17) 

Ohmic  heat  is 

rreff  ^-eff  Keff  Vct 

Qohm  =  -yr  •  Vfc  +  -^-V02  '  V02  +  -g--^  •  V4>2  (18) 

■ i  i  i  Z 

Active  heat  Qact  and  ohmic  Qohm  are  irreversible,  while  reaction 
heat  Qrea  is  reversible. 


Qirr  —  Qact  +  Qohm  (19) 

Ore  =  Qrea  (20) 


The  temperature-dependent  open  circuit  potential  of  electrode 
i  is  approximated  by  Taylor’s  first  order  expansion  around  a  refer¬ 
ence  temperature: 

Q;  =  Qref,i  +  (T-Tref)^  (21) 

where  Lfref,z  is  the  open  circuit  potential  under  the  reference  tem¬ 
perature. 

According  to  Newton’s  cooling  law,  the  boundary  condition  for 
energy  balance  is  expressed  as 

=  h(T,mb-T)--X^  =-h(7amb-r)  (22) 

x=0  x=Ln  +Ls+Lp 

where  h  is  lumped  heat  transfer  coefficient  and  Tamb  is  the  ambient 
temperature. 


3.3.  Solution  method 


The  calculations  of  model  were  performed  using  the  finite 
element  commercial  software  COMSOL  MULTIPHYSICS®  (Version 
3.5a).  Battery  performance  was  calculated  by  given  battery  cur¬ 
rent.  The  outputs  of  the  model  are  the  cell  potential,  current 
density  distribution,  species  and  concentrations  distributions.  The 
cell  potential  is  derived  by  the  following  expression: 

ECell  —  01 1  —  0i  I  (23) 

cen  ^ 1  \x=Ln+Ls+Lp  ^  1  lx=0 


4AA.  Model  parameters 

Model  parameters  are  from  experimental  measurement,  liter¬ 
ature  and  estimation.  Parameters  for  an  11. 5 Ah  LiMn204  battery 
are  listed  in  Table  1,  including  design  specifications,  lithium  ion 
concentration  parameters,  kinetic  and  transport  properties  and 
thermal  properties. 

4.  t.2.  Liquid  phase  parameters 

The  data  of  lithium  ion  diffusivity  in  solution  phase  D2,  electric 
conductivity  in  solution  phase  /c2  and  thermodynamic  factor  v  are 
adopted  from  the  published  literature  [9]: 

D2  —  lx  io-410“4-43(54-0/(r“229-0_0-05c))“2-2xlo“4c2  (24) 


k2  =  lx  10“4C2  (  -  10.5  +  0.074T  -  6.69  x  l(T5r2 

+  6.68  x  10“4c2  -  1.78  x  10“5c2T  +  2.8  x  10“8c2T2 
+  4.94  x  10 ”7cf  -  8.86  x  10"10c|t)2  (25) 


v  =  0.601  -0.24Vl0-3c2 


+  0.982[1 


0.0052(T  -  294.0) 


(26) 


Table  1 

Model  parameter  for  an  11.5Ah  LiMn204  HEV  battery. 


Anode 

Separator 

Cathode 

Design  specifications  (geometry  and  volume  fractions) 

Acdi  (m2) 

0.42753 

£u 

0.5a 

- 

0.49a 

£2,i 

0.33a 

0.54a 

0.332a 

Li  dun) 

120a 

30a 

150a 

Rj  (|Jim) 

12. 5a 

- 

8.5a 

Lithium  ion  concentrations 

C2,o  (molm-3) 

1200a 

Qnax,i  (tnol  m  3 ) 

26390b 

- 

22  860b 

SOCo,, 

0.74c 

- 

0.35c 

Kinetic  and  transport  properties 

<*a ,i,  ac,i 

0.5 

- 

0.5 

Yi 

1.5 

4d 

1.5 

D2  (m2  s_1) 

Eq.  (24) 

D10>I  (m2  s-1) 

3.9  x  10-14b 

- 

7.51  x  10“14d 

Ea,D,i  (kj  mol1 ) 

35e 

- 

31.556d 

k0j  (m2-5  mol-0-5  s-1) 

1.764  x!0-llf 

3.626  x  10-lld 

Ea,k.i  (kj  mol-1 ) 

20c 

- 

32.694d 

/ci  (Sm^1) 

100b 

- 

3.8b 

k2  (Sm4) 

Eq.  (25) 

t+ 

0.363b 

V 

Eq.  (26) 

F(Cmol4) 

96487 

Thermal  properties 

ki  (W(mK)-1) 

1.04s 

Is 

1.48s 

Pi  (kgm-3) 

2500h 

1200h 

1500h 

CpiO(kgK)-1) 

700h 

700h 

700h 

Tref(/0 

298.15 

a  Battery  manufacturer. 
b  Ref.  [9]. 
c  Ref.  [10]. 
d  Estimated. 
e  Ref.  [11]. 
f  Ref.  [6]. 
s  Ref.[12]. 
h  Ref.  [13]. 


Cell  potential  (V) 
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4 A. 3.  Kinetic  and  transport  properties 

Both  diffusivity  coefficient  and  reaction  rate  constants  are  fitted 
into  the  following  Arrhenius  formula  [10]: 


where  Di  n  and  /<i  n  data  for  MCMB  negative  electrode  come  from 
the  literature  as  shown  in  Table  1.  Battery  pulse  test  performance, 


Fig.  4.  Estimated  parameters  Dp  and  kp  are  fitted  into  Arrihenius  formula. 


as  shown  in  Fig.  2,  varies  at  different  temperatures,  mainly  due  to 
the  variation  of  Dlp  and  /<lp  at  different  temperatures.  The  reason 
may  be  discussed  later.  Regarding  to  this  property,  Dip  and  /<i  p 
are  estimated  and  fitted  in  Arrhenius  formula  by  comparing  model 
simulated  data  with  experimental  data.  The  results  are  presented 
in  Figs.  3  and  4. 

4 A. 4.  Lumped  heat  transfer  coefficient 

Lumped  heat  transfer  coefficient  h  is  determined  by  compar¬ 
ing  experimental  data  with  simulated  result  (Fig.  5)  .  Lithium  ion 
battery  is  discharged  with  2C  constant  current  with  temperature 
rising  to  50  °C  at  the  end  of  discharge.  Then  the  current  is  cut  off 
for  temperature  relaxation  to  room  temperature.  By  comparing 


0  2000  4000  6000  0  2000  4000  6000 


Time  (s)  Time  (s) 


Fig.  3.  Estimated  parameters  Dp  and  kp  by  comparing  1 C  pulse  test  experiment  data  with  simulated  data. 
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Time  (s) 


Fig.  5.  Determination  of  lumped  heat  transfer  coefficient  h  by  comparing  experi¬ 
mental  data  with  simulated  results. 


y  in  Li_Mn204 


0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 


Fig.  6.  Open  circuit  voltage  of  negative  electrode  Un  and  positive  electrode  Up. 


y  in  LLMn204 

0.0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 


x  in  Li_Cc 

x  6 


Fig.  7.  Entropy  change  of  negative  electrode  dUn/dr  and  positive  electrode  dUp/dT 
as  a  function  of  state  of  charge. 


Fig.  8.  Comparison  of  simulated  results  with  experimental  data  at  constant  currents 
(0.2C,  0.5C,  1C,  2C)  at  constant  temperature  25  °C. 


Discharge  capacity  (Ah) 

Fig.  9.  Comparison  of  simulated  results  of  battery  surface  temperature  with  exprei- 
mental  data  during  quasi  insulating  discharge  (0.2C,  1C,  2C). 


Charge  capacity  (Ah) 

Fig.  10.  Comparison  of  simulated  results  of  battery  surface  temperature  with 
expreimental  data  during  quasi  insulating  charge  (0.2C,  1C,  2C). 
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Fig.  11.  Comparison  of  measured  surface  temperature  between  quasi-insulating 
charge  and  discharge. 

temperature  relaxation  curve,  we  determine  the  heat  transfer  coef¬ 
ficient  to  be  0.38  W(m-2  K_1).  Lumped  heat  transfer  coefficient 
here  represents  the  synthetically  effect  of  heat  conduction,  con¬ 
vection  heat  transfer  and  radiative  heat  transfer. 

4.1.5.  Electrode  thermodynamic  properties 

As  shown  in  Fig.  6,  open  circuit  voltage  (OCV)  of  positive  elec¬ 
trode  (Up)  is  experimentally  (LAND  battery  test  system,  China) 
measured  using  method  present  in  Ref.  [11];  open  circuit  voltage 


(OCV)  of  negative  electrode  (Un)  comes  from  Ref.  [12].  In  Fig.  7, 
entropy  change  in  positive  electrode  is  A5P  =  nF(dDp/dT),  dUp/dT 
is  calculated  from  data  of  Ref.  [13]  using  method  presented  in 
Ref.  [14];  negative  electrode  entropy  change  is  ASn  =  nF(dDn/dF), 
dUnldT  data  comes  from  Ref.  [6]. 

State  of  charge  of  both  negative  and  positive  electrodes  has  the 
same  expression  as  follows: 

SOQ  =  Cu  (29) 

^1,  max,  i 


Fig.  13.  0.2C  pulse  test  validation  by  comparing  experiment  data  with  simulated  results. 
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That  is  SOCn  =  x  =  ciin/ci)max,n  in  LixC6  and  SOCp  =y  =  CiiP/ci)max,p 
in  LiyMn204.  Data  of  Up,  Un,  dUn/dT  and  dUp/dT  are  fitted  as  func¬ 
tions  of  SOC  by  piecewise  cubic  method  in  FUNCTION  dialog  box  of 
COMSOL. 

4.2.  Model  validation 

Model  validation  is  conducted  by  comparing  simulated  results 
with  experiment  data  at  different  conditions.  This  model  is 
validated  from  three  aspects:  electrochemically,  thermally  and 
electrochemical-thermally. 

4.2.  t.  Discharge  behavior  at  constant  temperature  condition 

Fig.  8  shows  the  comparison  between  simulated  results  and 
experimental  data  during  discharge  under  the  operating  condition 
at  constant  current  discharge  (0.2C,  0.5C,  1C,  2C)  at  constant  tem¬ 
perature  (25  °C).  The  results  indicated  that  the  modeling  results 
agreed  well  with  the  experimental  results  with  the  maximum  error 
less  than  2%. 

4.2.2.  Heat  sources  and  surface  temperature 

Fig.  9  shows  comparison  of  simulated  results  against  exper¬ 
imental  data  of  battery  surface  temperature  during  discharging 
process  in  quasi  insulating  condition.  It  can  be  seen  that  the 
modeling  results  of  battery  surface  temperature  at  both  0.2C  and 
2C  discharging  agreed  well  with  the  experimental  results,  while 
there  is  a  little  larger  deviation  at  the  end  of  1C  discharge  curve. 
The  results  indicate  that  the  battery  temperature  rise  faster  with 
larger  discharging  current.  When  the  discharge  current  is  large, 
the  temperature  increase  fast.  For  example,  when  discharging  at 


2C  discharge  current,  the  temperature  increase  about  30  °C  from 
the  beginning  at  discharging  process  to  the  end  of  the  discharging 
process.  In  this  situation,  the  battery  cooling  system  must  be  care¬ 
fully  designed  to  avoid  the  battery  safety  problem.  It  also  can  be 
seen  obviously  that  the  battery  surface  temperature  rise  faster  at 
the  beginning  of  discharging  process,  while  rise  slower  when  the 
discharge  capacity  is  larger  than  about  6Ah  (half  of  the  full  bat¬ 
tery  capacity).  When  at  the  beginning  of  the  discharging  process, 
the  temperature  is  relatively  lower,  the  battery  polarization  is  rel¬ 
atively  larger,  thus  the  more  heat  generated  during  this  period.  In 
addition,  the  lumped  heat  transfer  is  less  during  the  beginning  of 
discharging  since  the  temperature  difference  between  the  battery 
surface  and  environment  is  relatively  small. 

Fig.  1 0  also  shows  good  agreement  of  temperature  between  sim¬ 
ulation  and  experiment  during  charging  (0.2C,  1 C,  2C).  We  find  that 
surface  temperature  change  during  charging  and  discharging  are 
almost  the  same,  see  Fig.  1 1 ,  which  means  that  the  amount  of  heat 
generated  in  battery  during  low  rate  (0.2C,  1 C,  2C)  charging  and  dis¬ 
charging  are  almost  the  same  value.  We  will  further  discuss  heat 
generation  during  higher  rate  charging  and  discharging  below. 

4.2.3.  Dynamic  performance 

Fig.  1 2  shows  experimental  data  for  pulse  test  with  1 C  current  at 
different  temperatures.  During  one  cycle  of  discharge,  the  battery 
is  first  discharged  with  1 0  A  constant  current  for  360  s  and  then  cut 
off  the  current  for  relaxation  for  1800  s.  The  battery  is  discharged 
with  this  repeating  process  until  cell  potential  drops  to  2.5  V.  Dur¬ 
ing  the  discharge  period,  there  is  significant  cell  potential  drop, 
there  are  several  reasons  account  for  this  phenomenon.  First,  OCV 
of  battery  decrease  in  the  discharge  process.  Then,  over  potential, 
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Fig.  15.  2C  pulse  test  validation  by  comparing  experiment  data  with  simulated  results. 


consisted  of  concentration  polarization,  ohmic  potential  drop  and 
active  polarization,  accounts  for  the  increasing  gap  between  cell 
potential  and  OCV  of  battery.  At  point  A,  we  can  see  that  cell  poten¬ 
tial  drop  increases  with  the  temperature  drop,  due  to  lower  lithium 
ion  diffusivity  in  both  solid  and  liquid  phases,  larger  resistance  in 
battery  and  lower  reaction  rate  at  low  temperature.  During  relax¬ 
ation  period,  A-B,  active  polarization  and  ohmic  potential  drop 
extinct  and  result  in  sharp  increase  of  cell  potential.  Then  the  con¬ 
centration  polarization  takes  charge  and  leads  to  slower  increase  of 
cell  potential  until  B  spot.  After  1 800  s  relaxation,  cell  potential  at  B 
spot  is  approximately  the  same  as  OCV  of  the  battery.  Cell  potentials 
at  point  B  increase  as  temperature  increases,  and  this  phenomenon 
results  from  the  entropy  change,  see  Eq.  (21).  It  should  be  noticed 
that  entropy  change  is  not  always  the  same  at  different  SOC  (state 
of  charge)  of  battery,  so  the  OCV  regulation  may  not  necessarily 
applicable  for  other  SOC. 

Parameters  related  to  cell  dynamic  response  and  pulse  relax¬ 
ation  behavior  (e.g.  Dlp,  klp)  are  mainly  estimated  by  comparing 
simulated  results  with  experimental  data  at  1C  discharge  at  differ¬ 
ent  temperatures  (0°C,  10°C,  25  °C,  55  °C),  the  simulated  results 
for  pulse  test  (0.2C,  0.5C,  2C)  at  different  temperatures  (0  °C,  1 0  °C, 
25  °C,  55  °C)  were  validated  by  experimental  results  to  further 
check  the  model  applicability. 

Figs.  13-15  show  the  validation  results  at  0.2C,  0.5C  and  2C 
rates  at  different  temperatures  (0  °C,  1 0  °C,  25  °C,  50  °C).  The  results 
shown  that  the  modeling  results  agreed  well  with  the  experimental 
results,  while  have  a  little  larger  deviations  at  the  end  of  pulse-test 
at  each  temperatures.  The  deviations  are  preliminary  deduced  to  be 
due  to  the  changes  of  the  SEI  resistance  and  contact  resistance  with 
temperature  or  battery  state  of  charge  [15,16]  which  were  assumed 
as  constants  (zero)  in  our  model. 


5.  Results  and  discussion 

5.1.  Dynamic  response  of  lithium  ion  concentration  in  solid  phase 

One  of  the  potential  application  fields  of  the  proposed  model  is 
to  predict  the  lithium  ion  concentration  distributions  within  bat¬ 
tery  electrode  and  electrolyte,  which  is  helpful  for  battery  design 
and  performance  optimizations.  As  an  example,  pulse  test  with 
2C  discharge  rate  was  chosen  to  show  the  lithium  ion  concen¬ 
tration  distribution  response.  During  experimental  operation,  the 
battery  is  discharged  at  2C  constant  current  during  0-180  s  and 
then  is  relaxed  by  cutting  off  current  during  180-1000  s.  One  elec¬ 
trode  particle  at  boundary  1  was  chosen  to  examine  the  lithium 
ion  concentration  dynamic  change  at  different  positions  in  radius 
direction  (Fig.  16).  At  0-1 80s,  the  lithium  ion  concentration  gra¬ 
dient  increases  with  discharge  time,  lithium  ion  concentration  at 
r/Rn  =  0,  0.2,  0.4  decrease  during  0-1000  s,  while  lithium  ion  con¬ 
centration  at  r/Rn  =  0.6,  0.8, 1  decrease  during  0-1 80  s  and  increase 
during  180-1000  s.  Large  concentration  gradient  in  solid  particles 
may  lead  to  decrease  of  battery  capacity  and  battery  power  dur¬ 
ing  discharging,  which  will  be  discussed  below.  Reducing  particle 
radius  and  increasing  lithium  ion  diffusivity  in  solid  phase  would 
increase  the  relaxation  performance  of  battery  regarding  to  the 
relaxation  time  factor  r2 ID,  which  is  introduced  in  Ref.  [11]. 

5.2.  Effect  of  temperature  on  lithium  ion  distribution 

Fig.  17  shows  lithium  ion  concentration  distribution  across  the 
cell  in  liquid  phase  in  20  A  (2C)  pulse  tests  (t=180s).  During  dis¬ 
charge  process,  lithium  ion  concentration  in  negative  electrode  is 
higher  than  initial  concentration  and  lithium  ion  concentration  in 
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180s 
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Fig.  16.  Lithium  ion  concentration  distribution  along  particle  radius  during  pulse 
tests. 


Fig.  18.  Lithium  ion  concentration  in  solid  phase  along  particle  radius  at  different 
operating  temperatures. 


positive  electrode  is  lower  than  initial  concentration.  Since  lithium 
ion  is  extracted  from  negative  electrode  and  diffuses  to  liquid  phase, 
lithium  ions  intercalate  to  positive  electrode  particles  and  lower 
the  liquid  phase  lithium  ion  concentration.  The  lithium  ion  concen¬ 
tration  gradient  increases  with  decrease  in  operating  temperature. 
This  is  mainly  because  that  lithium  ion  diffusivity  increases  with 
the  operating  temperature. 

Fig.  18  shows  the  solid  phase  lithium  ion  concentration  dis¬ 
tribution  within  a  particle  at  current  collector/anode  interface. 
The  lithium  ion  concentration  gradient  in  particle  radius  direction 
increases  as  the  operating  temperature  decrease. 

After  discussing  the  lithium  ion  concentration  gradient  in  both 
solid  and  liquid  phases  upon  different  temperatures,  we  further 
investigate  the  battery  performance  variation  due  to  tempera¬ 
ture  change  (Fig.  19).  In  the  electro-thermal  model,  the  discharge 
current  is  set  to  20  A,  and  battery  is  discharged  at  different  temper¬ 
atures.  Battery  capacity  (product  of  discharge  time  and  discharge 
current)  decreases  significantly  with  temperature  drop,  mainly  due 
to  large  concentration  gradient  in  both  solid  and  liquid  phases  as 
discussed  above.  Battery  power  (integral  of  IxVxt)  also  decreases 
sharply  with  temperature  drop.  So  it  is  necessary  to  develop 
less  temperature-dependent  electrolyte  and  electrode  materials  to 
improve  battery  performance  at  lower  temperature,  since  high  per¬ 
formance  of  battery  at  low  temperature  is  required  for  winter  use 
in  the  application  of  EV/FIEV. 


5.3.  Reversible  heat  generation  and  irreversible  heat  generation 

In  this  section,  the  effect  of  entropy  change  during 
charge-discharge  on  battery  temperature  change  was  stud¬ 
ied.  As  mentioned  above  in  Eqs.  (19)  and  (20),  heat  generation  can 
be  divided  into  reversible  heat  Qre  and  irreversible  heat  Qirr.  Influ¬ 
ence  of  current  applied  to  battery  on  Qre  and  Qirr  is  approximated 
as  follows: 

Qirr  _  Qact  +  Qohm  ^  ^2^cell  ^  r  ^cell  (301 

Qre  “  Qrea  ~  IT(dU/dT)  ~  T(dU/dT )  ^  ; 

where  T  does  not  change  significantly  during  discharging  process, 
dU/dJ  is  a  function  of  state  of  charge  but  do  not  change  with  current 
/,  and  Rce\\  is  supposed  to  be  constant,  so  Qre  and  Qirr  are  dominated 
by  dU/dJ  and  Kcell,  respectively.  As  shown  in  Eq.  (30),  Qirr/Qre  is  pro¬ 
portion  to  /.  Fig.  18  depicts  different  parts  of  heat  generated  during 
charge-discharge  at  0.2C  rate.  Irreversible  heat  Qirr  and  reversible 
heat  Qirr  are  of  the  same  order.  Qirr  is  approximately  constant  with 
discharge  capacity,  while  Qirr  changes  with  discharge  capacity  sig¬ 
nificantly.  As  a  result,  the  battery  surface  temperature  curves  are 
of  the  same  shape  of  reversible  heat  Qre  in  charge  and  discharge 
processes,  respectively,  see  Figs.  20  and  21. 

Figs.  22  and  23  present  different  heat  sources  (reaction  heat  Qrea, 
ohmic  heat  Qohm*  active  polarization  heat  Qact)  during  0.2C  and  2C 
rate  discharge.  Irreversible  heat  (Qhom  and  Qact)  in  2C  discharge 


Negative  electrode  Separator  Positive  electrode 


Distance  from  negative  current  collector  (pm) 


Fig.  17.  Lithium  ion  concentration  distribution  within  liquid  phase  across  the  cell 
at  different  operation  temperatures. 


Fig.  19.  Influence  of  temperature  on  battery  2C  discharge  performance. 
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Fig.  20.  Different  parts  of  heat  generation  in  lithium  ion  battery  during  0.2C 
charge-discharge. 


Discharge  capacity  (Ah) 

Fig.  21.  Simulated  and  experimental  data  of  battery  surface  temperature  at  0.2C 
charge-discharge  rate. 

is  approximately  100  times  that  of  0.2C  discharge,  and  reversible 
heat  (Qrea)  in  2C  discharge  is  approximately  10  times  that  of  0.2C 
discharge,  this  phenomenon  fit  Eq.  (30)  quite  well.  Reversible  heat 
and  irreversible  heat  are  of  the  same  order  during  0.2C  discharge, 
while  irreversible  heat  is  about  10  times  that  of  reversible  heat 
during  2C  discharge,  based  on  this  phenomenon  and  Eq.  (30),  we 


Discharge  capacity  (Ah) 

Fig.  22.  Different  heat  sources  during  0.2C  discharge. 


can  safely  conclude  that  reversible  heat  dominant  at  low  rate  dis¬ 
charge  and  irreversible  heat  dominant  at  high  rate  discharge,  and 
we  can  ignore  reversible  heat  at  high  rate  (such  as  IOC)  discharge 
for  simplification. 

5.4.  Simulation  of  temperature  during  high  rat  charge  and 
discharge 

In  Fig.  24,  simulation  results  of  high  rate  (5C,  7C,  IOC) 
charge-discharge  are  shown.  We  consider  the  battery  is 
charged/discharged  with  proper  cooling  condition  with  a  bound¬ 
ary  condition  of  h  =  10Wm_2K_1,  each  temperature  curve  rises 
at  the  beginning  of  charging/discharging,  and  then  reaches  equi¬ 
librium,  the  same  as  Fig.  9  shows.  Different  from  Fig.  11,  Fig.  24 
shows  that  equilibrium  temperature  during  charging  is  higher  than 
that  of  during  discharging,  indicating  that  there  would  be  more 
heat  generated  during  charging  than  discharging  at  higher  C  rate. 
It  should  also  be  noted  that  temperature  rises  up  to  70  °C  at  the 
end  of  IOC  charging/discharging,  this  temperature  would  be  harm¬ 
ful  for  battery  and  could  be  dangerous.  More  intensive  cooling 
measure  should  be  taken  for  cooling  the  battery  during  high  rate 
charge-discharge. 


Fig.  24.  Simulation  results  on  temperature  change  during  high  rate  charge  and 
discharge. 
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6.  Conclusions 

In  this  study,  a  mathematical  model  by  coupling  electronic 
conduction,  mass  transfer,  energy  balance  and  electrochemical 
mechanism  is  developed.  Lithium  ion  diffusivity  and  chemical 
reaction  rate  of  cathode  material  are  estimated  by  comparing  sim¬ 
ulated  results  with  experimental  data  of  pulse  test  at  various 
current  charge-discharge  rates  (0.2C,  0.5C,  1C,  2C)  and  operating 
temperatures  (0  °C,  1 0  °C,  25  °C,  55  °C).  The  modeling  results  are  fur¬ 
ther  validated  in  aspects  of  electrochemical  performance,  thermal 
performance  and  electrochemical-thermal  coupling  effects,  which 
show  good  agreement  between  the  modeling  results  and  experi¬ 
mental  results. 

The  modeling  results  indicated  that  the  lithium  ion  concen¬ 
tration  gradient  in  both  liquid  phase  and  solid  phase  are  greatly 
affected  by  temperature  variation.  The  results  also  indicated  that 
the  lithium  ion  concentration  gradient  increases  when  the  operat¬ 
ing  temperature  decreases.  The  above  phenomenon  results  in  the 
capacity  losses  and  power  losses  of  lithium  ion  battery  during  low 
temperature  operation.  Then,  the  influence  of  reversible  heat  on 
battery  discharge  is  evaluated.  Reversible  heat  is  found  to  be  dom¬ 
inant  heat  source  during  low  rate  discharge  and  is  considered  to 


be  negligible  for  model  simplification  during  high  rate  discharge. 
Temperature  change  during  high  rate  charging/discharging  would 
be  too  high  for  battery,  proper  cooling  methodologies  should  be 
taken. 
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